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1 .  INTRODUCTION  AND  SUMMARY 


This  report  describes  how  to  implement  into  an  existing  satellite  lines 
status  computer  code,  the  author's  current  state  of  understanding  of  the  structuring 
and  convection  of  ionization  from  a  high-altitude-nuclear-explosion  and  the 
resultant  scintillation  of  radio  signals  penetrating  the  structure.  The  existing 
code  must  be  or  use  R ANC  IV  Mod  3 ,  a  nuclear  effects  code  built  by  G.  E. 

Tempo  for  the  D.  N.  A.  * 

The  changes  are  simple  and  straightforward  and  can  easily  be  implemented 
in  other  codes  which  use  the  RANC  phenomenology  routines. 

These  codes  currently  use  a  description  of  ionization  structure  obtained 

from  the  radiance  profile  of  a  single  late  time  nuclear  burst.  This  structure  is 

time -independent  and  has  a  power  spectrum  that  decreases  exponentially  in  wave- 

number.  It  is  currently  believed  that  the  power  spectrum  of  the  ionization  of 

typical  bursts  will  have  a  power  law  form  which  leads  to  greater  scintillation 

2 

effects  at  SHF  and  above  than  would  be  due  to  an  exponential  form. 

It  must  be  stressed  that  the  changes  that  we  are  making  in  the  structure 

description  are  not  the  result  of  a  complete  and  consistent  treatment  of  the  physics 

of  the  structuring  process.  That  treatment  has  not  yet  been  done  by  anyone  in  the 

community  of  contractors  to  the  DNA.  What  has  been  done  by  this  community  is 

a  numerical  nonlinear  analysis  of  the  structuring  of  ionization  with  a  one-dimen- 

3  4 

sional  Gaussian  variation  in  the  direction  of  a  constant  neutral  wind.  ’ 

We  explain  in  Section  2  how  the  results  of  this  numerical  analysis  lead 
to  insights  into  the  structuring  to  be  expected  in  the  ionization  of  a  late -time 
nuclear  burst. 

From  these  insights  we  suggest  a  model  of  the  time  dependent  structuring 
of  the  ionization  which  specifies  parameters  of  the  structure  that  relate  directly  to 
scintillation  descriptions.  These  parameters  do  not  uniquely  follow  from  the 


3 


numerical  analysis.  There  may  be  alternate  models  and  parameters  that  are  a 
better  representation  of  the  structure.  We  view  the  modeling  described  herein  as 
a  first  cut  to  demonstrate  the  feasibility  of  the  concept  of  this  type  of  approach  to 
the  fast-running  scintillation  model  and  invite  other  members  of  the  community  of 
contractors  to  devote  some  of  their  time  to  its  improvement. 

When  structuring  begins,  the  ionization  is  able  to  convect  across  magnetic 
field  lines.  Current  codes  do  not  include  this  effect  but  restrain  the  ionization  by 
magnetic  containment.  We  use  the  results  of  the  same  numerical  analysis  to  model 
this  convection. 

In  Section  3  we  describe  how  we  implement  the  structure  and  convection 

g 

model  into  the  RANC  fireball  model  using  the  results  of  a  hydro  code  calculation 
to  guide  the  correspondence  between  the  numerical  analysis  results  and  the  nuclear 
case.  In  Section  4  we  describe  the  scintillation  effects  in  terms  of  the  parameters 
of  structure  and  ionization  content  using  extensions  of  results  obtained  in  a  previous 
study.  ® 

In  Section  5  we  describe  the  suggested  method  of  implementing  the  struc¬ 
ture  and  scintillation  model  into  the  existing  code.  A  by-product  is  a  change  in  the 
calculated  absorption  to  be  consistent  with  the  changed  ionization  model. 

The  essential  concepts  of  the  stucture  and  convection  model  are:  There 
is  a  period  of  time  after  a  high  altitude  burst  when  expansion  across  the  ambient 
magnetic  field  has  ceased  and  the  ionization  is  constrained  by  the  magnetic  field. 

We  use  the  time  TEQ  from  the  RANC  code  for  this  period. 

After  this  time  neutrals  continues  to  have  a  radial  component  of  velocity 
across  the  magnetic  field  causing  the  gradient  drift  instability  to  be  operative. 

The  results  of  a  numerical  calculation  of  the  effects  of  that  instability  on  the 
structuring  and  convection  of  a  one -dimensionally  varying  plasma  is  applicable 
for  scaling  the  initial  structuring  and  convection  of  the  nuclear  plasma.  The 
motion  of  the  plasma  in  the  direction  of  the  magnetic  field  and  the  time  variation 
of  the  electron  density  are  uncoupled  from  the  structuring  and  convection  across 
the  field  and  are  taken  from  the  RANC  code.  There  is  a  scaled  time  after  TEQ 
when  structuring  has  begun  internally  but  the  radius  to  the  boundary  of  the 
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ionization  as  given  by  RANC,  RTF,  has  not  changed.  After  this  time  the  boundary 
radius  grows  with  time  at  a  constant  rate  as  large  fingers  of  ionization  protrude 
radially  whose  tips  bifurcate  into  smaller  fingers  of  lesser  electron  density  until 
all  the  ionization  is  in  the  two  sets  of  fingers.  After  this  time  the  ionization  in 
the  large  fingers  flows  into  the  smaller  fingers  which  have  bifurcated  again. 
Finally,  the  small  fingers  continuously  bifurcate  until  a  minimum  size  is  reached. 
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2.  THE  STRUCTURE  MODEL 

The  structure  model  is  based  on  the  following  concepts.  First,  the  various 
two-fluid  hydrocodes  (e.  g. ,  MELT,  built  by  MRC)  are  useful  for  the  prediction  of 
the  early  time  (<  500  sec.)  distribution  of  ionization  following  a  high-altitude  nuclear 
burst  and  the  conditions  which  cause  structuring  (neutral  wind,  gravity,  etc.). 

They  cannot,  however,  describe  the  structuring  that  does  occur  because  of  the 
very  large  (~  100  km)  cell  size  they  must  use.  After  structuring  should  have 
occurred,  these  codes  will  be  incorrect  not  only  by  the  absence  of  the  structuring 
but  also  by  the  absence  of  the  convection  of  the  ionization  across  magnetic  field 
lines  because  of  the  structuring. 

Second,  an  illustration  of  the  structuring  and  convection  unobtainable 
from  the  hydrocodes  is  given  by  the  output  of  a  numerical  electrostatic  code 
such  as  one  existing  at  NRL. 

The  only  example  of  such  a  calculation  that  we  feel  is  pertinent  to  this 

nuclear  structuring  problem  was  done  for  the  case  of  a  plasma  in  a  uniform  mag- 
3  4 

netic  field.  ’  The  plasma  had  a  Gaussian  variation  of  density  in  one  direction 

superimposed  on  a  constant  background.  The  problem  was  done  for  the  case  of 

a  uniform  electric  field  at  right  angles  to  the  magnetic  field  and  to  the  density 

gradient.  This  is  completely  equivalent  to  the  case  of  a  uniform  neutral  wind 

blowing  through  the  plasma  at  right  angles  to  the  electric  and  magnetic  fields 

in  the  gradient  direction.  In  the  absence  of  perturbations  the  plasma  would 

remain  stationary  with  the  drag  of  the  neutral  wind  balanced  by  a  magnetic  force 

(3  x  §).  The  situation  at  the  beginning  of  the  calculation  is  illustrated  in  Figure  1 

3 

taken  from  the  NRL  report.  The  neutral  wind  is  in  the  y  direction  and  the 
current  that  balances  the  drag  is  in  the  x  direction  and  varies  only  with  y. 

The  pertinence  of  this  calculation  to  the  nuclear  case  is  the  high  ratio  of 
the  maximum  density  to  the  background  and  the  lack  of  variation  of  the  plasma  in 
the  direction  of  f.  In  the  nuclear  case  the  plasma  geometry  is  essentially  cylin¬ 
drical  with  the  axis  along  a  relatively  uniform  magnetic  field.  The  density 

i-riLCLui^G  f>.GE  blauiK-MOT  F IxJCG 
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Fig.  16.  Isodensity  contours  of 
plasma  density  at  t  —  0  sec.  The  ini¬ 
tial  destribution  for  Nf/N0  is  a  gaus- 
sian  in  y,  centered  at  y  —  12.1  km, 
plus  a  small  random  perturbation  in  x. 
Contours  are  drawn  for  NJN„  —  1.5, 
3.5,  5.5,  7.5  and  9.5.  The  area 
between  every  other  contour  line  is 
cross-hatched.  Only  120  of  the  160 
cells  actually  used  in  the  y  direction 
are  displayed.  Boundary  conditions 
are  periodic  in  both  directions.  In  our 
plot  B„  is  toward  the  reader,  and  E„  is 
directed  toward  the  right,  and  we 
have  placed  outselves  in  a  frame 
moving  with  the  (c/J^I2)  E0x  B„ 
velocity.  The  upper  portion  of  the 
gaussian  is  physically  unstable  to  per¬ 
turbations,  while  the  lower  half  is 
(linearly)  stable. 


OSEC 


X  (KM) 


Figure  1.  Initial  Contours  NRL  3716. 
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variation  and  the  neutral  wind  direction  is  radial.  The  current  is  azimuthal  and 
varies  in  the  radial  direction.  The  structuring  begins  near  the  edge  of  the 
plasma  where  a  cartesian  approximation  to  radial -azimuthal  variation  is  best. 

The  subsequent  behavior  of  the  edge  of  the  plasma  in  response  to  the  small  per¬ 
turbations  in  the  numerical  calculation  is  similar  to  the  behavior  expected  at  the 
edge  of  the  nuclear  plasma  during  the  first  phase  of  its  structuring.  This  is  illus¬ 
trated  in  Figure  2  taken  from  the  same  report.  These  figures  indicate  certain 
concepts  of  structure  that  are  useful  for  the  modeling  philosophy. 

Before  describing  these  concepts  we  will  review  the  insights  obtained 
from  analysis  of  the  simple  problem  of  an  elliptic  cylindrical  plasma  cloud  of 
constant  conductivity  with  axis  along  a  constant  magnetic  field  B  =  Bz,  im¬ 
bedded  in  a  background  plasma  of  conductivity  E^.  A  neutral  wind  of  uniform 
velocity  V  =  Vy  blows  through  the  plasma.  If  the  background  plasma  is  not  mov¬ 
ing  in  the  region  far  from  the  cloud,  there  is  no  electric  field  in  that  region.  The 
cloud  boundary  contains  a  polarization  surface  charge  producing  an  electric  field 
in  and  near  the  cloud.  The  motion  of  the  cloud  or  background  plasma  is  governed 
by  the  equation 


V  = 
P 


-*  -* 
E  X  B 


The  electric  field  is  found  from  the  equations 


v.  J  =  V  x  E  =  0 


J  =  E(E  +  V  x  B)  =  E(E  +  VBx)  . 

where  far  from  the  cloud,  Jo  =  E^VBx  .  There  is  a  discontinuity  in  the  normal 
component  of  E  at  the  cloud  boundary  because  of  the  discontinuity  in  £.  From  the 
divergence  theorem  J  •  n  =  J  •  n  or  £  (E  +  VBx)  •  n  =  £  (E  +  VBx)  •  n, 

0  D  C  C  Dt)  ^ 

where  n  is  a  unit  normal  at  the  boundary.  Adding  and  subtracting  •  n  obtains 
the  formula  for  the  surface  charge 
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(KM) 


r  *. 


tr  ■  (Efe  -  Ec)  •  n  =  <  ^!3)  J  •  n  (1) 

choosing  n  to  point  to  the  background. 

The  solution  for  the  electric  field  everywhere  is  obtained  by  the  usual 
procedure  for  a  dielectric  in  a  field  which  is  uniform  at  infinity.  Set  E  =  -  v<p 

o 

where  ?  <p  =  0  within  and  without  the  plasma  cloud  and  <p  -*  0  at  infinity.  At  the 
boundary,  <p  is  continuous  and  v<p  has  a  jump  according  to  the  surface  charge 
density  a. 

For  the  special  case  of  an  elliptic  cloud  with  ratio  of  y  intercept  to 
x  intercept  j3,  an  analytic  solution  shows  that  the  electric  field  within  the  cloud 
is  uniform  with  value 


->  -+ .  2 

The  cloud  moves  as  a  whole  with  velocity  Ec  x  B/B  .  For  the  usual  case  where 
f&c  »  2^,  the  velocity  is  Vc  **  Vy(l  -  At  the  boundary  E^*  n  =  -(1  n 

and  the  interior  current 

?  J.d  +  0  -  /  l 

The  insights  gained  from  the  elliptical  problem  are 

•  The  higher  the  conductivity  ratio,  the  more  the  cloud  moves 
with  the  neutral  wind. 

•  The  higher  the  elongation  in  the  wind  direction  the  less  the 
exterior  electric  field  and  change  in  current. 

•  If  0  »  1  and  2^  >  Ec,  the  cloud  moves  against  the  neutral 
wind  with  a  speed  greater  than  neutral  wind  speed  by  a 
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factor  j  |3/[1  +  tfE/iy ]  j . 

We  use  these  insights  to  describe  some  of  the  qualitative  behavior  shown  by  the 
conductivity  contours  in  Figure  2.  The  initial  configuration  was 

2 

£  =  1  +  10e“(y//8)  [1  +  e(x,  y)] 

-4 

where  c  was  generated  from  a  K  power  spectrum  with  random  phases  and  *rms 
=  0. 03.  There  is  a  neutral  wind  of  100  m/s  in  the  y  direction.  The  contours  are 
£  =  1.  5,  3.  5,  5.  5,  7.  5,  9.  5. 

From  the  two  figures  it  is  apparent  that 

•  Fingers  of  high  conductivity  move  slower  than  fingers  of  low  conduc¬ 
tivity  that  go  against  the  wind.  The  plasma  therefore  separates  into 
many  large  blobs  in  a  short  time. 

•  The  steepening  occurs  at  the  lower  conductivity  contours  in  the 
plasma  fingers.  This  is  because  the  change  in  electric  field  and 
hence  the  change  in  velocity  is  proportional  to  o  which  by  expression 
(1)  is  proportional  to  (£  -2y/£c  since  J  is  relatively  uniform. 

The  Gaussian  profile  in  y  becomes  broader  until  £~  5.  5  where  it 
steepens  abruptly.  The  x  profile  of  course  is  the  finger  width. 

•  Two  background  fingers  descended  with  little  resistance.  The 

high  conductivity  interior  and  the  lower  region  of  decreasing  con¬ 
ductivity  do  not  begin  to  move  until  the  fingers  pass  them.  Because 
of  the  periodic  boundary  conditions,  the  voltage  drop  across  the 
mesh  remains  zero,  i.  e. ,  <p(x  =  0)  =  =  L).  Therefore  the  lower 

regions  do  not  experience  an  electric  field  until  the  passage  of  the 
low  conducitivity  finger  with  its  own  electric  field  requires  an  ex¬ 
terior  electric  field  to  keep 

L 

/E  dx  =  0  . 
x 

O 
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The  behavior  of  the  third  background  finger  in  the  central  region  is  different 
since  it  is  not  in  a  uniform  background  of  high  conductivity.  Its  bifurcation 
produced  a  steep  small  high  conductivity  finger.  It  is  interesting  that  the 
other  plasma  fingers  appear  to  have  median  values  of  £  where  they  tend  to  form 
necks  prior  to  break  off.  This  signifies  small  striations  with  low  values  of 
central  conductivity  occurring  in  many  cases  with  the  bulk  of  the  peak  values 
of  conductivity  in  the  central  region  of  large  striations. 

By  about  450  sec. ,  the  NRL  run  ends  showing  two  main  highly  elongated 
fingers  of  width  3  km  and  two  smaller  fingers  of  smaller  length  and  width,  one 
high  density  and  one  at  lower  average  density.  The  trifurcation  that  appears  in 
the  one  finger  in  the  figure  also  begins  in  the  other  main  figure  at  450  sec.  We 
do  not  understand  the  tendency  to  trifuracate  in  the  NRL  calculation.  We  can 
speculate  that  the  tendency  of  the  finger  to  flatten  and  steepen  sets  the  stage  for 
the  next  phase  of  fingering  and  in  the  case  of  the  NRL  calculation,  a  0. 3  km  cell 
size  and  a  3  km  finger  width  limits  the  number  of  new  fingers  to  three.  The 
tendency  to  flatten  and  steepen  is  apparent  in  Figure  2.  The  finger  on  the  right 
has  moved  an  average  velocity  of  60  m/s  between  the  two  times  whereas  the 
£  =  5.  5  and  7.  5  contours  have  moved  an  average  of  80  m/s.  The  finger  on  the 
left  undergoing  trifurcation  has  moved  at  50  m/s.  The  boundary  of  the  plasma 
has  moved  12  km  during  the  calculation  time  for  an  average  initial  convection 
speed  during  the  first  phase  of  30  m/s. 

We  must  now  scale  the  results  of  the  NRL  calculation  to  what  we  expect 

to  occur  in  the  first  phase  of  structuring  in  the  ionization  from  a  high  altitude 

nuclear  burst.  The  first  question  we  address  is  the  size  of  the  fingers  to  be 

expected.  The  NRL  run  had  fingers  of  width  3  km.  Although  the  Gaussian  scale 

length  for  the  radial  £  dependence  was  8  km,  the  arguments  put  forth  in  Chapter 

6 

4  of  a  previous  analysis,  indicate  the  real  scale  length  in  question  is  the  mini¬ 
mum  value  of 

h  e  ivin£|  * 

which  is  4.  7  km  for  the  NRL  case.  A  3  km  finger  corresponds  to  a  wavelength  of 
6  km  so  that  defining  k  =  2tr/X,  we  have  kh^^  »  5  which  we  take  as  the  scaling 
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required  to  estimate  finger  width  in  the  nuclear  case. 

The  next  question  we  address  is  the  scaled  time  required  to  finish  the 
first  phase  of  fingering  which  was  about  450  sec.  in  the  NRL  case.  According 
to  the  linear  instability  theory,  the  growth  rate  =  1 A  *  V/hmin  or  Tg  "  ^  sec. 
Since  they  started  with  3%  perturbations,  the  nonlinear  phase  should  begin  at  t 
=  r  An  33  =  3.  5  r  =  165  sec.  This  occurs  in  the  narrow  region  of  maximum 

g  g 

growth  (minimum  h)  illustrated  in  Figure  3.  It  takes  a  little  longer  for  the  last 
contour  (2  =  1.  5)  to  resemble  the  contour  of  the  budding  fingers.  We  take  as  the 
time  for  the  background  to  have  fingered  in  sufficiently  to  begin  rapid  movement 

through  the  cloud  as  4  r  .  Now  the  problem  should  begin  to  resemble  the  previ- 

S 

ously  described  analytic  case  of  an  ellipse  in  a  constant  background  except  the 
background  is  not  constant.  In  that  case  an  ellipse  made  up  of  background  con¬ 
ductivity  flowing  through  a  constant  cloud  of  higher  conductivity  went  against  the 
neutral  wind  with  a  speed 


The  nonuniformity  of  the  background  should  reduce  this  speed  by  an  un¬ 
known  amount.  Referring  back  to  Figure  2  we  see  qualitative  agreement  in  that 
background  fingers  go  faster  than  plasma  fingers  but  the  two  background  fingers 
seem  to  be  travelling  at  ~  90  m/s  between  300  and  400  sec.  Comparison  of 
Figures  2  and  3  show  the  speeds  to  be  90  and  135  m/s  between  200  and  300  sec. 
Thus  the  nonuniformity  of  the  background  reduces  the  analytic  speed  expected 
considerably. 

Since  the  NRL  calculation  with  its  constraint  that  the  polarization  elec¬ 
tric  field  average  to  zero  over  the  x  dimension  of  the  mesh  should  be  more 
similar  to  the  nuclear  case  (where  the  azimuthal  electric  field  should  have  zero 
average)  than  the  analytic  case  of  uniform  infinite  bakcground,  we  choose  our 
scaling  of  the  background  finger  speed  from  the  NRL  case.  We  take  the  back¬ 
ground  finger  speed  as  the  neutral  wind  speed.  Thus  the  time  for  this  phase  will 
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108  SEC 


204  SEC 


Figure  3.  Intermediate  Contours  NRL  3716. 


be  the  time  for  the  fingering  to  begin,  4  r  ,  plus  the  time  for  a  background  finger 
to  break  through  at  the  neutral  wind  speed.  For  the  NRL  calculation  this  is  188 
sec.+  23.  5  km/0. 1  km/sec.  =  423  sec. 

During  the  235  sec.  after  fingering  has  begun,  the  plasma  finger  moved 
out  at  an  average  speed  of  half  the  neutral  speed,  50  m/s.  We  will  therefore  take 
the  convective  displacement  of  the  plasma  boundary  as  1/2  V(t  -  4  f^). 

What  happens  after  this  phase  ?  Another  NRL  calculation  of  the  same 

problem  was  carried  out  to  720  sec.  This  is  illustrated  in  the  R  =  “plot  in 

4  e 

Figure  4  taken  from  NRL  3864.  The  other  plots  for  finite  Rg  have  diffusion  added 
which  is  too  large  for  the  parameter  region  of  interest  here. 

We  find  three  main  fingers  with  an  average  convective  displacement  in 
agreement  with  our  formula.  The  boundary  at  the  right  is  at  50  km  and  one  finger 
has  passed  and  re-entered  on  the  left  (periodic  boundary  conditions).  A  second 

4 

has  reached  the  boundary  and  a  third  has  not.  This  calculation  solves  the  same 

3 

differential  equations  with  the  same  initial  conditions  as  the  other.  Slight  differ¬ 
ences  in  numerical  techniques  cause  slight  differences  in  the  contours  which 
are  close  enough  for  our  modeling  purposes.  The  interesting  phenomenon  in 
this  later  time  picture  is  the  persistence  of  the  original  three  fingers.  Their 
tips  have  become  small  fingers  of  lesser  conductivity  but  the  regions  of  high 
conductivity  remain  intact  in  the  large  fingers. 

It  has  been  said  that  after  this  phase  the  radial  fingers  in  the  nuclear 
case  will  be  in  a  uniform  neutral  wind  and  will  act  like  Ba  clouds.  Ba  clouds 
start  out  spherical,  flatten,  and  finger.  Only  the  tips  of  these  phase  1  fingers 
act  like  Ba  clouds.  We  have  no  numerical  calculations  or  Ba  experiments  that 
tell  us  what  happens  to  the  bulk  of  the  highly  elongated  fingers.  We  therefore 
conjecture  that  after  this  phase  the  boundary  of  the  ionization  continues  to  move 
out  at  half  the  neutral  wind  speed  and  is  made  up  of  two  sets  of  forms.  Each  of 
these  forms  will  be  described  using  the  four  parameter  description  of  a  stria- 
tion  whose  scintillation  effects  were  described  in  an  earlier  paper.  The  electron 
density  profile  is  taken  as 


* 


ttwait  jam** 


Figure  5 

Z  contours  at  T  ■  720  seconds 


Figure  4.  Later  Contours  NHL  3864 
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r  <  a 


=  0  r  >  a 

2  2  2 

where  r  =  0x  +  y  /0.  Contours  of  electron  density  are  ellipses  with  major  to 
minor  axes  ratio  0.  This  profile  is  parabolic  from  nQ  down  to  nj  =  nQ(  1  -a2) 
where  the  fall-off  is  abrupt.  The  four  parameters  nQ,  a,  %  and  0  describe  the 
striation. 

The  first  set  consists  of  the  large  striations.  The  peak  electron 
density,  nQ,  is  taken  as  the  peak  of  the  original  unstructured  ionization.  The 
semi-minor  axis  is  half  the  finger  width 

\  _  ^min 

The  specification  of  a  and  0  requires  a  consistency  in  preserving  the  con- 
o 

tent  of  n  and  n  in  addition  to  approximating  the  elongation  and  profile  of  the 
major  portion  of  the  content. 

The  same  is  true  for  the  second  set  which  represents  the  secondary  fingers. 
We  assign  the  electron  density  of  the  second  set  the  value  n^  =  nQ(l-c?)  which 
makes  it  dependent  on  the  first  set.  We  assign  the  value  a  =  0  for  the  second 
set  for  convenience.  The  initial  size  of  the  secondary  fingers  are  half  the 
primary.  The  semi-minor  axis  of  the  second  set  is  then 

as  _  X  _  ^min 

■irt  ’ 8 '  20  ' 

The  specification  of  a  and  the  elongation  factors  0T  and  0  set  the 

J-i  s 

content  in  each  striation.  The  large  striations  each  have  content 

"0*4  (x  ■  ®2/2) 
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V 


where  a£  =  0L(ffhmin/lO)2. 


The  small  strialions  each  have  content 


n^tra 


2 

s 


where  a2  =  M«h  .  /20)2  . 
s  s  min 

The  original  content  is  divided  between  the  two  requiring 


where  M_  and  M  are  the  number  of  large  and  small  striations,  a  is  the  original 
Li  S  § 

Gaussian  scale  of  the  unstructured  ionization  and  the  factor  '1  -  or)  in  the  second 

term  is  n1/n  .  The  original  unstructured  ionization  had  a  squared  ionization  of 

2  2  10  2  2  4 

n  ira  /2  of  which  division  between  the  two  requires  M  ira  (1  -  a  +  or /3)  +  M 

O  On  a  »  ri  LiLi  S 


2x2  2  2 

(1  -  a  )  ira  =  ira  /2  .  The  solution  of  the  equation  is 
s  o 


2 

which  requires  a  *  1/2.  This  arises  from  the  constraint  of  fixing  the  square  of 

n  content  of  a  Gaussian  by  means  of  parabolic  and  rod  profiles.  The  model 

e  2 
accomplishes  the  requirements  of  preserving  n  content  and  n  content  while 

modeling  the  large  sizes’  maintenance  of  smooth  variation  and  the  small  sizes 

existing  through  pinch-off  of  large  size  edge  bifurcation.  The  scintillation  effects 

of  the  small  sizes  vanish  at  a  =  1. 


In  order  to  relate  the  first  set  with  the  scintillation  effects  described  in 
2  2  2 

Reference  6  we  redefine  aT  ( 1  -  a  /2)  =  a  so  that  the  content  of  a  large  striation 
2  L 

is  nQiia  .  The  results  of  that  study  are  described  in  terms  of  a  rather  than  a^ 
The  original  content  n  ira^  is  thus  split  into  CT  and  C  where 

O  O  S 


2 

CT  =  n  M.  ira  =  n  ite 
L  o  L  oo 


2  (2a2  -  1)  (1  -  a2/2) 
a2(l  -  a2/3) 
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■  noraofs(o) 


C  =  M  ira2  (l  -  a2)n  =n 
s  s  s  o  i 


.  wa2 
o  o 


37  - 14 


Once  the  two  contents  are  established,  the  further  changes  in  striation 

2 

morphology  are  modeled  by  decreasing  a  and  aa  while  keeping  N,  a  and 

2  ®  "  2 

N  a  constant.  In  addition  further  pinch  off  is  modeled  by  decreasing  a 
s  s 

toward  l/2  thereby  increasing  C  at  the  expense  of  CT  . 

S  Ju 
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3.  IMPLEMENTATION  OF  STRUCTURE  AND 
CONVECTION  INTO  FIREBALL  MODEL 

The  structure  model  discussed  in  the  last  section  redistributes  the  elec¬ 
tron  content  per  length  along  the  magnetic  field  that  is  predicted  by  a  hydrocode 
such  as  MELT.  In  this  section  we  present  an  example  of  that  hydrocode  output 
and  illustrate  the  redistribution.  We  describe  the  output  as  a  structured  RANC 
fireball. 

Since  we  are  using  the  RANC  fireball  model  as  the  basic  descriptor  of 
the  unstructured  ionization  in  general,  we  will  use  the  MELT  output  as  an  inter¬ 
face  between  the  RANC  fireball  and  the  structure  model  since  MELT  supplies  us 
with  an  ionization  profile  whereas  RANC  does  not.  The  RANC  ionization  is  uni¬ 
form  in  the  plane  perpendicular  to  the  magnetic  field. 

MRC  has  kindly  supplied  us  with  curves  of  electron  density  from  a 
MICE/MELT  calculation  of  a  particular  high  altitude  nuclear  burst  which  is 
representative  of  yield  and  burst  height  to  be  expected.  Figure  5  is  an  example. 
We  have  fit  their  output  at  two  calculation  times,  500  and  1000  sec.  by  specify¬ 
ing  the  maximum  electron  density  nQ(z,  t)  (where  z  is  the  distance  along  the 

magnetic  field)  and  the  total  electron  content  per  unit  length  along  z,  C(z,  t). 

2  2 

The  content  is  n  ira.  where  a  is  the  product  of  two  characteristic  radii  from 
oo  o 

Gaussian  fits  to  the  electron  density  profiles  in  two  directions  perpendicular 
to  the  magnetic  field.  An  example  of  the  contours  of  electron  density  in  a 
plane  perpendicular  to  the  magnetic  field  is  shown  in  Figure  6.  The  MELT 
runs  show  an  asymmetry  which  we  ignore  in  using  the  symmetric  RANC  fire¬ 
ball.  We  do  not  believe  the  changes  that  would  be  required  in  the  RANC  inter¬ 
section  routines  to  accommodate  this  asymmetry  are  warranted. 

The  RANC  fireball  model  is  a  truncated  ellipsoid  described  by  the 
parameters: 


HF;  Height  of  center  (km) 


600  SEC 


NE 


<CM-3> 


CYCLE  2 

PHJ  r  -4  3€E-| 


Figure  5.  Contours  in  Plane  of  Magnetic  Field. 


(SNViova-miH) 
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RLF;  Semi-major  axis  (along  magnetic  field)  (km) 

RTF;  Semi-minor  axis  (km) 

TILTF;  Cosine  of  magnetic  dip  angle 
HMINF;  Height  of  fireball  bottom  (km) 

HMAXF;  Height  of  fireball  top  (km) 

Within  this  ellipsoid  RANC  prescribes  the  electron  density  as  a  function 
of  height  only. 

The  values  of  these  parameters  for  the  burst  and  times  corresponding  to 
the  two  MELT  runs  are: 


Time 

(sec) 

HF 

RLF 

RTF 

TILTF 

HMINF 

HMAXF 

500 

1695 

3199 

473 

0. 2474 

335.4 

4893 

1000 

1695 

6063 

473 

0.  2674 

22a  0 

7758 

The  resulting  shapes  are  sufficiently  close  to  the  MELT  contours  to  use 
them  to  contain  the  ionization  content  of  the  MELT  runs. 

The  agreement  is  good  when  the  RANC  ellipsoid  is  adjusted  to  lie  along 

the  Earth’s  dipole  magnetic  field  as  is  currently  done  in  G.  E.  Tempo's  WEPH 
7  8 

code  and  SAI's  code  NSM.  We  therefore  recommend  this  adjustment  in  conjunc¬ 
tion  with  our  structure  model. 

We  therefore  take  the  RANC  ellipsoid  as  the  boundary  of  the  ionization 
obtained  from  MELT.  We  will  be  identifying  this  boundary  with  the  results  of  the 
NRL  numerical  calculation  discussed  in  Section  2  by  associating  it  with  the 
boundary  that  was  above  the  background  and  illustrated  the  fingers.  The  RANC 
boundary  corresponds  to  electron  density  contours  two  orders  of  magnitude  down 
from  the  peak  MELT  values  but  higher  than  the  background. 

The  right  way  to  associate  the  NRL  calculation  with  the  MELT  output  is 
through  a  magnetic  field-line  integrated  Pederson  conductivity  and  a  conductivity- 
weighted  field-line  integrated  neutral  wind  and  other  forcing  terms  obtained  from 
the  MELT  output.  Since  we  do  not  have  these  quantities  we  cannot  accurately 
evaluate  the  radius  transverse  to  the  field  lines  where  the  log-conductivity  gradi¬ 
ent  maximizes.  We  therefore  choose  the  RANC  radius  as  qualitatively  reasonable. 
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Comparison  of  a  roughly  field -line -integrated  value  of  a  with  the  RANC 

O  A 

radius  at  500  sec.  shows  aQ  *  RTF/2.  From  a  previous  analysis,  we  know  that 

a  Gaussian  profile  of  conductivity  with  a  maximum  value  two  orders  of  magnitude 

over  the  background  has  a  minimum  gradient  length  =  aQ/3  which  is  only 

slightly  different  for  higher  maximum  to  background  ratios.  We  therefore  take 

h  .  =  RTF/6, 
mm 


Typical  values  of  the  radial  neutral  wind  speed  in  hydrocode  output  run 
from  100  m/s  at  a  transverse  radius  of  100  km  to  1000  m/s  at  the  fireball  edge.  ® 
In  Section  2  we  estimated  a  time  of  four  growth  times  for  the  structuring  to  be¬ 
come  evident  where  the  growth  time,  r  =  h  .  /V  where  V  was  uniform.  For 

g  mm 

the  nuclear  case  V  varies  linearly  with  radius  and  is  about  750  m/s  in  the  region 
of  maximum  initial  growth. 


In  the  RANC  fireball  model,  the  early  radial  expansion  of  the  plasma 
stops  at  a  time,  TEQ.  At  this  time  after  burst  when  they  fix  the  fireball  radius 
RTF,  the  plasma  is  at  time  zero  for  the  NRL  numerical  calculation.  We  there¬ 
fore  set  the  time  for  structuring  to  be  observable  as 


T 


1 


=  TEQ  +  4  t  =  TEQ  + 
g 

=  TEQ  +  (Sec. ) 


a  RTF(km) 

3  0.  75  km/s 


For  t  <  we  consider  the  fireball  to  be  unstructured. 

For  t  >  Tj,  the  first  phase  of  structuring  proceeds  until  the  background  fingers 

have  penetrated  to  the  center  of  the  fireball.  In  Section  2  this  time  was  determined 

from  a  constant  velocity.  For  the  nuclear  case  of  radial  penetration  of  background 

fingers,  the  velocity  has  a  linear  variation  with  radius  which  leads  to  a  logarithmic 

infinite  time  of  penetration.  Since  the  correspondence  of  the  numerical  calculation 

with  the  nuclear  case  is  best  at  the  plasma  edge,  we  make  the  analogy  of  finger 

displacement  at  the  plasma  edge  where  the  plasma  fingers  (and  hence  the  boundary) 

convected  out  an  amount  D  =  \  V(t  -  r^).  The  size  of  the  fingers  is  determined 

from  the  scaling  law  kh  .  =  5  or  k  =  30/RTF  since  h  .  =  RTF/6.  As  the  plasma 

mm  mm 

fingers  move  out,  the  background  must  move  in  to  preserve  the  area  in  the  original 
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circle  of  radius  RTF.  The  background  finger  will  have  the  form  of  a  wedge  (pie 

o 

slice)  of  area  (RTF)  X/4  =  ir(RTF)  /60  when  the  background  finger  reaches  the 
origin.  At  this  time  the  plasma  finger  will  have  an  area  beyond  the  original  circle 
corresponding  to  a  half  ellipse  or  ffD/2  X/4  which  will  equal  the  pie  slice  when 
D  =  2  RTF/ if  at  the  time 


r 


2 


,  4  RTF 
ffV 


where  V  **  1  km/s. 


At  the  end  of  this  time  the  boundary  of  the  striated  ionization  will  be  at 
RTF  +  2  RFT/ it.  The  structure  is  described  by  the  parameters  mentioned  in 
Section  2.  The  parameter  a  distributes  the  content  C  between  the  two  sets  of 
sizes,  determines  the  profile  of  the  large  sizes,  and  the  relative  electron 
density,  n/n  ,of  the  small  sizes.  The  large  sizes  have  maximum  electron 

O  t 

density  nQ(z).  We  shall  begin  at  t  =  t  with  a  =  1  and  end  at  t 2  with  n^/nQ  =  (1-a 
=  1/20  which  is  the  ratio  of  conductivity  where  the  gradient  length  h  is  smallest 
to  the  maximum  conductivity.  At  the  same  time  we  will  be  turning  on  the  struc¬ 
ture  with  a  coefficient  that  linearly  rises  from  zero  at  to  unity  at  to  repre¬ 
sent  the  structuring. 


Ct  Z  Ct 

For  a  "  1  we  have  from  Section  2,  CT  **  t  n  ra  and  C  "»in  ifa  . 

L  o  o  s  o  o 

For  the  initial  number  of  large  fingers  we  use  the  criterion  k  =  30/RTF 
which  signifies  the  number  of  fingers,  M^,  to  be  30  since  that  is  the  number  of 
wavelengths  around  the  circumference  of  the  ionization.  Therefore  equating  the 
formulas  for  content,  30  nna~  /2  =  f  nna  or  aT  =  a  /20  =  (RTF)  /80.  From 

n  O  Li  0  0  a  ij  O  a 

Section  2  a,  =  0T  (ffh  .  /10)2  =  0T  (irRTF/60)  .  Therefore  0T  =  (60/n)  /80  ~  4. 6. 
i-j  Lt  min  L  l 


Since  the  small  set  represents  the  bifurcation  of  the  tips  of  the  large 

fingers,  their  number  and  elongation  are  not  constant.  Also,  since  we  are  com- 

2 

mited  to  conserving  the  n  and  n  content  of  a  Gaussian  profile  using  non- 

0  0 

Gaussian  profiles  we  suffer  the  consequence  that  the  initial  content  of  the  small 

fingers  be  one  quarter  the  total  which  is  not  true.  Their  initial  scintillation 

2 

effects  however  are*n  and  start  from  zero.  We  shall  set  the  initial  size  of 

e 


26 


the  small  fingers  so  that  their  width  is  half  that  of  the  large  fingers  so 


We  expect  /3g  **  0^  =  4. 6  so  we  will  let  their  number  be  determined  from  the  content 
formulas.  Choosing  the  large  size  as  aL  =  RTF/9  says  the  initial  small  size  is 
ag  =  RTF/18. 

For  t  >  Tg,  we  expect  the  striations  to  continue  their  radial  convection 
but  further  structuring,  if  any,  will  not  be  apparent  until  after  the  last  time  of 
the  numerical  calculation  shown  in  Figure  4  of  Section  2 .  Because  of  the  high 
elongation  (fl)  and  high  conductivity  ratio  to  background  and  lower  slip  velocity 
(neutral  speed-cloud  speed)  we  expect  further  structuring  to  take  a  relatively 
long  time.  It  would  seem  reasonable  that  the  smaller  fingers  continue  to  grow  at 
the  expense  of  the  large  fingers.  The  larger  fingers'  higher  conductivity  will 
flow  into  the  small  fingers. 

By  allowing  a  to  approach  %  we  will  deplete  the  large  fingers.  We 
expect  further  bifurcation  to  occur  in  the  small  fingers.  The  typical  way  this 
occurs  need  not  result  in  a  much  smaller  finger  width  since  the  fingers  can  widen 
from  the  filling  before  bifurcating.  Thus  as  a  -»  4,  the  electron  density  in  the 
small  fingers  and  their  number  increase  and  their  size  decreases  by  only  half. 

We  have  no  numerical  calculations  or  experimental  data  that  are  late  enough 
to  fix  the  time  for  the  end  of  this  phase  when  a  =  i.  We  shall  take  this  time 
as  the  time  for  the  high  density  (n  =  nQ/2)  to  overtake  the  boundary.  The  high 
density  is  taken  as  i  the  original  peak  since  this  is  the  final  value  of  density  in 
the  small  finger  when  a  =  By  the  analysis  in  Section  2,  the  high  density 
portion  should  travel  at  the  local  neutral  wind  speed  while  the  boundary  will 
continue  to  travel  at  half  of  the  local  neutral  wind  speed.  At  t  =  t^,  the  boundary 
is  at  RB  =  RTF(1  +  2 /if)  and  the  high  density  is  at  RH  =  0.  6  RTF.  We  take  the 
neutral  wind  speed  as  V  =  MIN  (r/RTF,  1)  km/s  which  is  the  speed  of  the  high 
conductivity  portion.  The  high  density  position  a  time  t  later  is 

RH  =  RTF  +  t  -  RTF  In  (5/3)  =  0.  49  RTF  +  t 
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j£  * 


where  the  last  term  in  the  middle  is  the  time  to  reach  r  =  RTF. 

RB  =  RTF  ^1  + 

=  1.637  RTF  +0.  5t  . 

We  find  RH  =  RB  at  t  =  2.  29  RTF  after  We  therefore  define 

t3  =  t2  +  2.  29  RTF 

when  we  are  left  with  only  the  small  set  of  striations  described  as  elliptical 

"waterbags"  (0  =  4. 6,  a  =  0)  of  ionization  n1  =  n  /2  and  size  a  =  RTF/36. 

1.  o  s 

If  they  were  truly  waterbags  they  would  not  structure  further  if  they 
were  in  a  uniform  neutral  wind  according  to  the  analysis  of  Section  2.  We  know 
that  any  departure  from  these  idealized  conditions  cause  them  to  structure  fur¬ 
ther  until  they  reach  some  stable  size  whose  value  has  not  yet  been  determined 

for  the  nuclear  case.  At  this  stage  the  striations  are  similar  to  barium  clouds. 

0 

In  a  previous  analysis  we  invoked  a  bifurcation  time  which  depended  on  the 
structure  parameters  and  conductivity  ratio  X-  Since  this  time  was  proportional 
to  the  striation  size,  each  bifurcation  time  was  progressively  shorter.  From 
barium  experiments  there  appeared  to  be  a  minimum  size  where  bifurcation 
ceased.  The  reason  for  this  has  not  yet  been  discovered  so  scaling  to  the 
nuclear  case  cannot  be  done  at  this  time.  Assuming  therefore  the  nuclear  case 
will  be  similar  to  the  barium  case  we  chose  am  =  200  m.  We  implement  this 
by  the  formula  from  reference  6  which  says  the  time  to  reach  size  a  from  an 
initial  size  a'  is 


2\  t 
i/  2 


km 


where 


t(a,  a')  =  r^i 


16  a 


a  -  a  + 


m 


a  -  a 


m 


16  a 


_  4. 3  «/ff 
°"  Vn 
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let 


and 


then 


t(a,a') 

T 

O 


16  a 


„  -  m 

1  a  -  a 


m 


a  = 


1 

2 


+  a 


m 


+  a 


m 


We  choose  the  parameters  x  =  4,  0  =  4.  6,  Vn  =  0.  5  km/s  and  get  tq  =  553  sAm. 

So  for  time  t  after  t 3  when  aQ  =  a'  -  (t  -  t3)/553,  am  =  0.  2,  and  aj 
-  0.  64/(a  -  0. 2)  where  a  =  RTF/36  we  use  the  above  formula  for  a 

s’ 

To  summarize: 


Given  RTF  (km)  and  TEQ  (sec. )  from  the  RANC  code; 

Define  Tj  =  TEQ  +  8/9  RTF  (sec.), 

T2  =  Ti  +  4/ffRTF  (sec.), 
and  t3  =  t2  +  2*  29  RTF  (sec. )  . 

Define  a  structure  factor  S  as  a  multiplier  for  the  scintillation  effects 
to  be  introduced  later.  The  time  independent  parameters  of  structure  are  0 

~  ~  ^  ~  RTF/9,  a^  =  0.  2  km.  The  time  dependent  structure  -ameters 

are  a  and  ag.  The  boundary  of  the  striations  is  at  RB. 

For  t  <  Tj 


S  =  0,  RB  =  RTF  . 


For  Tj  <  t  *  t2 


_  RTF 
s  18 


S  = 


t  -  r. 
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(t  -  T  ) 

,  RB  =  RTF  + - s-^- 


For  r2  <  t  *  t3 


as  = 


RTF 

TT 


1  - 


(t  -  T2) 


,  S  =  1  , 


2  1  9  -  T_)  (t  -  T  ) 

“  =1-20-^  (rpp  ’  =  RTF  +  g"” 


For  t  >  Tr 


S  =  1  , 


RB  =  RTF  + 


(t  -  Tj 


a  =  i  |a  -a  -  a„  +  /fa  -  a  )2  -  2a*  (a  +  a  -  2a  *)  +  a? ! 
s  I  o  m  1  \/v*o  m;  lv  o  m  }  1 ) 


where 


(‘  -  r3) 
553 


16  a 


„  m 

a< - r - 

1  a  -  a 


m 


,  _  RTF 
a  36 


The  final  portion  of  the  fireball  model  is  the  stipulation  of  the  unstruc- 
2  2 

tured  content  C  =  n  tra„  =  n  ff/4  (RTF)  and  the  peak  ionization  n  .  The  RANC 
ooo  o 

model  outputs  an  ionization,  ENE,  which  is  a  decreasing  function  of  distance 

2 

from  the  bottom  along  the  magnetic  field  but  is  uniform  within  the  area  HR 
where  R  is  the  transverse  radius  of  the  RANC  ellipsoid  with  semi-minor  axis 
RTF.  At  that  point 

ENE  w(RTF)2  =  nQw(RTF)2/4 
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or 


n  =  4  ENE  . 
o 


We  chose  aQ  =  RTF/2  as  an  approximation  for  a  field  line  integrated 
average.  By  using  the  2-D  NRL  run  we  are  saying  the  mixing  and  convection 
of  the  plasma  occur  by  rearrangement  of  the  contents  of  flux  tubes.  This 
means  that  while  the  plasma  boundary  is  fingering  where  R  =  RTF,  nothing  is 
happening  up  the  field  lines  where  there  is  only  background  at  R  =  RTF.  This 
means  the  higher  regions  should  be  retarded  in  their  structuring  compared  to 
the  lower  regions.  We  cannot  model  this  effect  now.  Since  we  must  work  with 
the  RANC  elliptic  fireball  we  simply  make  the  substitution  of  RB  for  RTF  and 
let  the  upper  part  of  the  ellipse  expand  accordingly.  We  maintain  the  parameters 
of  structure  throughout  the  fireball.  For  any  point  along  the  field  line,  the  RANC- 
predicted  unstructured  content  is 

2  /  z2 

ENE  fl(RTF)  1 - - — * 

\  (RLFj 

We  will  take  nQ  =  4  ENE  throughout  thereby  implying  that  the  Gaussian  scale  also 
decreases  with  z  in  the  same  manner.  The  implication  of  this  choice  is  that  the 
number  of  striations  per  unit  area  are  independent  of  z. 
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4.  SCINTILLATION  EFFECTS 


The  scintillation  effects  of  the  structure  model  were  determined  in 

Reference  6  for  a  single  size.  The  scintillation  effects  for  the  two  sets  are 

2 

obtained  by  adding  their  defined  in  that  reference  extended  to  the  case  where 
the  ray  path  is  not  perpendicular  to  the  magnetic  field  but  is  at  an  angle  6. 

D 


where 


„  16a3r.  rXn  2 

J  =  _  *i\  V  °'  N  P 

voL  3  sin  &  WL^ 


where 


r  is  a  scale  factor  caused  by  striation  ellipticity  which  will  be  described  later. 

r  =  2.  818  •  10"13  cm. 
e 

X  =  signal  wavelength. 

Nl=  number  of  striations  per  unit  area  perpendicular  to  the  magnetic  field. 

P  =  ray  path  length  through  the  striations. 

d  =  djdg/(dj  +  dg)  and  dj  and  dg  are  the  two  distances  between  the  two  terminals 
and  the  striations. 

The  striations  lie  within  an  area  w(RB)  in  the  plane  perpendicular  to 
the  magnetic  field.  Therefore 
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which  is  a  scale  factor  derived  in  Reference  6  where  0.  is  the  angle  between  the 
long  axis  of  the  elliptic  striation  and  the  projection  of  the  LOS  in  the  plane  per¬ 
pendicular  to  the  magnetic  field. 

As  a  typical  LOS  penetrates  the  striated  region  there  will  be  a  range  of 
values  of  6.  and  hence  r  will  be  a  weighted  average  which  will  be  near  V?  for  a 
LOS  crossing  near  the  center  of  the  striations  and  near  l/«/£  for  a  LOS  passing 
near  the  endge  of  the  striated  region.  As  explained  in  Reference  6,  this  occurs 
because  the  striations  have  their  long  axes  pointing  radially  from  the  center  of 
the  striated  region  in  the  plane  perpendicular  to  the  magnetic  field.  This  is 
modeled  by  choosing  8.  =  cos 


(P  sin  0g/2  RB). 
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The  correlation  length  of  the  scintillations 


2T^J\  -  o2/2 
o 

according  to  the  single  size  analysis  of  Reference  6.  The  extension  of  the  analysis 
to  two  sizes  defines  l  by  the  equation 


7 


=  4 


I^L(1  -  7/2) 


a 


a. 


2 

The  quantities  and  t  are  the  primary  descriptors  of  the  scintillation 
effects.  Scintillation  in  the  received  signal  intensity  is  measured  by  the  scintil¬ 
lation  index  defined  by 


-  1 


where  I  is  the  signal  intensity  at  the  receiver  normalized  to  intensity  in  the 
absence  of  structure  in  the  ionization.  The  average  intensity  T  =  1  by  conserva¬ 
tion  of  energy.  From  the  analysis  of  Reference  6, 


S2  -  1  -  e'2”e“ 
4 


2 

meaning  intensity  scintillation  is  negligible  for  <p&{f  «  1  and  essentially  Rayleigh 
for  >  1.  When  Rayleigh,  l  is  the  autocorrelation  length  of  I  -  F.  In  general 
l  is  the  length  such  that  E(x)E*(x  +  i)  =  lA/e  =  0. 606  where  E  =  ^1<P  where 
the  phase  <p  is  produced  by  the  structure  and  x  is  distance  perpendicular  to  the 
plane  defined  by  the  LOS  and  the  magnetic  field.  Therefore  in  the  absence  of 
intensity  scintillation,  l  has  meaning  for  phase  scintillation 


ei[tfx)  -  tfx  +  i)]  =  e-$[tfx)  -  < tfx  +  _  e-l 
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i  wwiiihi 

_ _  mm* 


jfc—.A.;..,.. 


or  fv<x)  -  <p(x  +  i)]2  =  1  . 

The  translation  of  correlation  length  to  correlation  time  requires  the 
x  component  of  velocity  of  the  LOS  with  respect  to  the  striations.  This  is 
sufficient  if  the  striations'  motion  is  negligible.  We  have  the  striations  moving 
across  the  magnetic  field  at  speeds  around  500  m/s.  If  the  LOS  passes  near 
center,  very  little  striation  motion  is  in  the  x  direction.  If  the  LOS  passes 
near  the  edge  of  the  striated  region  the  striation  motion  is  mainly  in  the  x 
direction.  We  therefore  take  V  =  0.  5  km/s  sin  6.  where  0.  was  used  in  the 

X  11 

definition  of  r  and  choose  the  sign  according  to  whether  the  LOS  is  moving 
radially  outward  or  inward.  That  is,  the  correlation  time 


T  = 


w. 


xLOS 


l _ 

±0.5  sin  0^  I 


This  formula  erroneously  allows  t  to  become  infinite.  We  know  the 
differential  movement  of  the  striations  will  cause  t  to  always  be  finite.  The 
determination  of  r  for  differentially  moving  striations  is  a  task  for  future 
work. 
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5.  IMPLEMENTATION  INTO  AN  EXISTING  CODE 


In  any  existing  systems  code  that  uses  RANC  phenomenology,  there  is  a 
routine  that  describes  the  fireball.  In  RANC  the  routine  is  called  PHENOM.  In 
that  routine  the  type  of  high-altitude  fireball  that  our  structure  model  is  pertinent 
to  is  magnetically  contained  so  that  after  a  defined  magnetic  containment  time 
TEQ,  the  semi-minor  axis  of  the  ellipsoidal  fireball,  RTF,  is  held  constant  at 
the  magnetic  containment  radius  RM.  In  PHENOM,  this  type  of  fireball  is 
designated  by  setting  the  flag  KINDF  =  2.  We  save  the  radius  RM  for  use  as 
RTF  in  our  structure  formulas  and  change  the  boundary  by  setting  RTF  =  RB 
according  to  the  formulas  of  Section  3  thereby  expanding  the  fireball  to 
account  for  the  striation  convection. 

There  is  also  a  routine  that  computes  the  attenuation  and  scintillation 

effects.  In  RANC,  the  routine  is  called  ABSORB.  It  calls  a  routine,  XFBT, 

that  determines  the  intersection  of  the  line  of  sight  between  two  antennas  and 

the  fireball.  It  then  steps  along  the  total  path  length,  PL,  of  the  intersection 

2 

adding  up  the  absorption  in  dB  and  the  mean-square  scintillating  phase, 
according  to  its  own  structure  formula.  These  two  quantities  are  proportional 
to  the  square  of  the  local  electron  density,  ENE,  which  it  obtains  from  the  sub¬ 
routine  ENEFB.  At  this  point  we  change  ENE  to  correspond  to  the  mean  square 
electron  density  of  our  expanded  fireball.  Therefore 


since  n  =  4  ENE. 
o 

For  the  absorption  formula  we  replace  ENE  by  JL  ENE(RM/RTF). 
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r*  » 


2  2 

For  scintillation  effects  we  need  nQ  =  16(ENE)  at  each  step  to  obtain 


o  _  64  Sa  -p  t 

'  T5-  L  L 


(reXRM\2 
\  RTF  / 


(ENE)*(  AP){ 
sin  0„ 


steps 


W.  ,„2 

°L 


The  determination  of  T  =  T  described  in  Section  3  requires  information 

S  Jj 

unobtainable  from  the  RANC  intersection  routine,  XFBT,  without  altering  the 
routine.  We  therefore  approximate  the  average  angle  between  the  LOS  and  the 
striation  axis  as 


0.  =  cos 

l 


\steps 


AP  sin  6 

O 

2  RTF 


for  the  formula 


r  =  rr  =1 
s  L 


1  +  (/3  -  1)  sin  0. 


The  formulas  for  a,  a  ,  P,  f  ,  fT  ,  ot,  and  S  given  in  Section  3  are  implemented 

S  S  l-i 

to  evaluate  and  output  and  i  from  the  formulas  in  Section  4  . 

The  RANC  code  outputs  the  velocity  of  the  LOS  perpendicular  to  the  mag¬ 
netic  field.  What  is  needed  is  the  component  of  that  velocity  in  the  x  direction 
and  the  position  of  the  LOS  intersection  in  the  fireball  in  order  to  compute  correla¬ 
tion  time.  Without  this  information  we  suggest  bounding  the  correlation  time 
using  V  =  500  m/s  added  to  and  subtracted  from  the  given  LOS  velocity  with  a 
minimum  value  of  50  m/s  for  the  denominator  in  the  formula  t  =  1/  IV^Qg  ±  0.  5 1  . 
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